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A new generalization of the P\ non-conforming FEM 
to higher polynomial degrees* 

M. Schedensackt 


Abstract 

This paper generalizes the non-conforming FEM of Crouzeix and Raviart 
and its fundamental projection property by a novel mixed formulation for 
the Poisson problem based on the Helmholtz decomposition. The new for¬ 
mulation allows for ansatz spaces of arbitrary polynomial degree and its 
discretization coincides with the mentioned non-conforming FEM for the 
lowest polynomial degree. The discretization directly approximates the gra¬ 
dient of the solution instead of the solution itself. Besides the a priori and 
rnedius analysis, this paper proves optimal convergence rates for an adaptive 
algorithm for the new discretization. These are also demonstrated in numer¬ 
ical experiments. Furthermore, this paper focuses on extensions of this new 
scheme to quadrilateral meshes, mixed FEMs, and three space dimensions. 

Keywords non-conforming FEM, Helmholtz decomposition, mixed FEM, adaptive FEM, 
optimality 

AMS subject classification 65N30, 65N12, 65N15 


1 Introduction 

Non-conforming finite element methods (FEMs) play an important role in com¬ 
putational mechanics. They allow the discretization of partial differential equa¬ 
tions (PDEs) for incompressible fluid flows, for almost incompressible materials in 
linear elasticity, and for low polynomial degrees in the ansatz spaces for higher- 
order problems. The projection property of the interpolation operator of the Pi 
non-conforming FEM, also named after Crouzeix and Raviart |21j . states that the 
L 2 projection of V-ffg (^) onto the space of piecewise constant functions equals the 
space of piecewise gradients of the non-conforming interpolation of Hq(£1) func¬ 
tions in the Pi non-conforming finite element space. This property is the basis for 
the proof of the discrete inf-sup condition for the Stokes equations [2T] as well as 
for the analysis of adaptive algorithms [5]. 

Many possible generalizations of the Pi non-conforming FEM to higher poly¬ 
nomial degrees have been proposed. All those generalizations are either based on 
a modification of the classical concept of degrees of freedom [23l|22] HD , are re¬ 
stricted to odd polynomial degrees EDI E|, or employ an enrichment by additional 
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bubble-functions |29l [28]. However, none of those generalizations possesses a cor¬ 
responding projection property of the interpolation operator for higher moments 
(see Remark [3TT5] below). This paper introduces a novel formulation of the Poisson 
equation (in (13.31) below) based on the Helmholtz decomposition along with its dis¬ 
cretization of arbitrary (globally fixed) polynomial degree. This new discretization 
approximates directly the gradient of the solution, which is often the quantity of 
interest, instead of the solution itself. For the lowest-order polynomial degree, the 
discrete Helmholtz decomposition of [3J proves equivalence of the novel discretiza¬ 
tion with the known non-conforming Crouzeix-Raviart FEM |21| and therefore 
they appear in a natural hierarchy. In the context of the novel (mixed) formu¬ 
lation, these discretizations turn out to be conforming. Although the complexity 
of the new discretization itself is competitive with that of a standard FEM, the 
method requires the pre-computation of some function ip such that its divergence 
equals the right-hand side. If this is not computable analytically, this results in an 
additional integration (see also Remark 13.61 belowh However, this paper focuses on 
the Poisson problem as a model problem to introduce the idea of the new approach 
and to give a broad impression over possible extensions as quadrilateral discretiza¬ 
tions (including a discrete Helmholtz decomposition on quadrilateral meshes for 
the non-conforming Rannacher-Turek FEM [33] as a further highlight of this pa¬ 
per), the generalization to three dimensions, or inhomogeneous mixed boundary 
conditions. The advantages of the new approach in some applications will be the 
topic of forthcoming papers mm- 

The presence of singularities for non-convex domains usually yields the same 
sub-optimal convergence rate for any polynomial degree. This motivates adaptive 
mesh-generation strategies, which recover the optimal convergence rates. This 
paper presents an adaptive algorithm and proves its optimal convergence. The 
proof essentially follows ideas from the context of the non-conforming Crouzeix- 
Raviart FEM j6j [32] . This illustrates that the novel discretization generalizes it 
in a natural way. Since the efficient and reliable error estimator involves a data 
approximation term without a multiplicative power of the mesh-size, the adaptive 
algorithm is based on separate marking. 

A possible drawback of the new FEMs is that the gradient of the solution Vu is 
approximated, but not the solution u itself. This excludes obvious generalizations 
to partial differential equations where u appears in lower-order terms. 

The remaining parts of this paper are organized as follows. Section [2] defines 
some notation. Section [3] introduces the novel formulation based on the Helmholtz 
decomposition and its discretization together with an a priori error estimate. The 
equivalence with the P\ non-conforming FEM for the lowest-order case is proved in 
Subsection 13.31 Section [4] summarizes some generalizations. Section [5] is devoted to 
a medius analysis of the FEM, which uses a posteriori techniques to derive a priori 
error estimates. Section [6] proves quasi-optimality of an adaptive algorithm, while 
Section [3 outlines the generalization to 3D. Section [8] concludes this paper with 
numerical experiments. 
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2 Notation 

Throughout this paper 0 C M 2 is a simply connected, bounded, polygonal Lip- 
schitz domain. Standard notation on Lebesgue and Sobolev spaces and their norms 
is employed with L 2 scalar product Given a Hilbert space X, let 

L 2 {Vt\X) resp. H k (fl]X) denote the space of functions with values in X whose 
components are in L 2 (Q) resp. H k (Fl) and let Tq(O) denote the subset of L 2 (Q) 
of functions with vanishing integral mean. The space of L 2 functions whose weak 
divergence exists and is in L 2 is denoted with iL(div, fl). The space of infinitely 
differentiable functions reads C' 00 (n) and the subspace of functions with compact 
support in hi is denoted with C%°(£1). The piecewise action of differential operators 
is denoted with a subscript nc. The formula A < B represents an inequality 
A < CB for some mesh-size independent, positive generic constant C\ A ~ B 
abbreviates A < B < A. By convention, all generic constants C ~ 1 do neither 
depend on the mesh-size nor on the level of a triangulation but may depend on 
the fixed coarse triangulation To and its interior angles. The Curl operator in two 
dimensions is defined by Curl/3 := (d/3/dx2, —dff/dx i) for sufficiently smooth /3. 

A shape-regular triangulation T of a bounded, polygonal, open Lipschitz do¬ 
main Q C M 2 is a set of closed triangles T £ T such that = (J T and any two 
distinct triangles are either disjoint or share exactly one common edge or one ver¬ 
tex. Let £(T) denote the edges of a triangle T and £ := £(T) := Up eT £(T) the 
set of edges in T. Any edge E £ £ is associated with a fixed orientation of the unit 
normal ve on E (and te = (0, —l;l,0)z /e denotes the unit tangent on E). On 
the boundary, ve is the outer unit normal of hi, while for interior edges E % dfl, 
the orientation is fixed through the choice of the triangles T + £ T and T_ € T 
with E = T + 0 T_ and ve '■= v t + \e is then the outer normal of T + on E. In this 
situation, [v\e '■= v\ t + — v\t_ denotes the jump across E. For an edge E C dQ on 
the boundary, the jump across E reads [v\e ■= v. For T £ T and X C M n , let 

each component of v is a polynomial 1 
of total degree < k J ’ 

P k ( T; X) := {v : Q -> X | VT £ T : v\ T £ P k (T ; X)} 

denote the set of piecewise polynomials and P k { T) := P k { T;M). Given a subspace 
X C L 2 (fi;M n ), let : L 2 (fI;M n ) —> X denote the L 2 projection onto X and let 
n k abbreviate Hp^y-j-igny Given a triangle T £ T, let Ht '■= (meas 2 (T)) 1 / 2 denote 
the square root of the area of T and let hj £ Pq(T) denote the piecewise constant 
mesh-size with hj \ t := hr for all T £ T. For a set of triangles M C T, let || • Hm 
abbreviate 


II • ll| 2 (T)- 

T£M 

Given an initial triangulation To, an admissible triangulation is a regular triangu¬ 
lation which can be created from To by newest-vertex bisection |40| . The set of 
admissible triangulations is denoted by T. 



Pk(T-X ) := {v.T^X 
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3 Problem formulation and discretization 

This section introduces the new formulation based on the Helmholtz decomposition 
in Subsection 13. II and its discretization in Subsection 13.21 Subsection 13.31 discusses 
the equivalence with the Pi non-conforming Crouzeix-Raviart FEM [21]. 


3.1 New mixed formulation of the Poisson problem 

Given the simply connected, bounded, polygonal Lipschitz domain fl C R 2 and 
/ € P 2 (fl), the Poisson model problem seeks u £ Hq(CL) with 

—A u = f in fl and u = 0 on dCt. (3-1) 


The novel weak formulation is based on the classical Helmholtz decomposition [35] 

P 2 (f!;R 2 ) = XHq(Q) © Curl(P' 1 (P) fl Lo(fl)) (3.2) 

for any simply connected domain fl C R 2 , where the sum is orthogonal with 
respect to the L 2 scalar product. 

Remark 3.1. Note that for fl C R 2 , the definition of the Curl implies 

H (Curl,fl) := {/3 £ L 2 (fl) | Curl/3 £ P 2 (fl)} = H l {Tt). 

Define X := T 2 (fl;R 2 ) and Y := nLo(fl) and let £ P(div,fl) satisfy 

— div</? = /. The novel weak formulation of the Poisson problem m seeks 
(p, a) £ X x Y with 


(P> q)l\S1) + (?> Curla) i2(Q) = (tp, q) L 2 (Q) for all q £ A, 

(p, Curl (3) l 2 (q) = 0 for all ft £ Y. 

This formulation is the point of departure for the numerical approximation of Xu 
in Subsection 13.21 


Remark 3.2 (existence of solutions). Since CurlT C X, any (5 £ Y satisfies the 
inf-sup condition 


ll Curl ^llL2( n) 


(g,Curl/3) L 2 (n) 

< sup -n n- 

gex\{o} \\q\\L2(n) 


This and Brezzi’s splitting lemma m imply the unique existence of a solution 
(p,ot) £ X x Y to (13.311 . The L 2 orthogonality of p and Curia implies 

\\p\\L 2 (Vt) + ||Curl a:||(j}) = WtWlz^)- 


Remark 3.3 (equivalence of (13.11) and (13,31) 1. The second equation of (13.31) and 
the Helmholtz decomposition (13.2p imply the existence ofu£ Hq (fl) with p = Xu. 
Since ip £ P(div,fl) satisfies — div (p = f, the L 2 orthogonality in (13.21) implies 
that any v £ Hq(£ 1) satisfies 


(P,Vu) L 2 (f2) = (t,Xv) l2(q) = (f,v) L 2(n) 
and, hence, u solves m- 
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Remark 3.4 (mixed boundary conditions). Let dQ = U Tjy with Tjj closed, 
r^nTiv = 0, and each connectivity component of To has positive length. Assume 
that the triangulation resolves To- Let H^ l ^ 2 (T^) denote the space of generalized 
normal traces of H(div,Tl) functions and let uo G H 1 (fi) and g G H~ 1 / 2 (Tn) 
in the sense that there holds g = q ■ v on Tn in the sense of distributions for 
some q G H( div, fi). Consider the mixed boundary value problem —A u = f in 
Tl with = U D on To and (Vu • i / )|r JV = g on T at. Let H^(Tl) denote the 
subspace of H 1 ^) of functions with vanishing trace on To- For To = 9, define 
H}j(n) := H 1 ^) n L 2 (Q). Define 

Hl(Q) := {/3 G Y | fd is constant on each connectivity component ofTjy}. 

The Helmholtz decomposition 

L 2 (Q]R 2 ) = ©Curlfli(n) 

for mixed boundary conditions \2f\ Corollary 3.1] then leads to the following for¬ 
mulation. Let <p G H( div, 0) with — div<£> = / additionally fulfil the boundary 
condition ipv\r N = g and seek ( p,a ) G L 2 (fi;R 2 ) x Hf(Tl) with 

(p,Q)L2(n) + (q, Curl a) L 2 ^ = (<p,q) L 2 ( n) for all q G L 2 (Q]R 2 ), 

(p, Curl/?)£ 2 (q) = (VuD,Curl/3) L 2 (n) for all /3 G 

Since p = tp— Curia G R(div, fi), the equivalence follows as in Subsection \3.1\ and 
with 


ip ' ^)|rv = (T ■ v)\r N - (Curia • v)\r N = g - (Va • t)\ Tn = g. 

Remark 3.5 (multiply connected domains). If Tl C R 2 is a multiply connected 
polygonal bounded Lipschitz domain and dTl = T/j U r^r, such that all parts of To 
lie on the outer boundary ofTl (on the unbounded connectivity component ofR 2 \Tl), 
then the Helmholtz decomposition of Remark \3.4\ still holds and a discretization as 
above is then immediate. However, if the Dirichlet boundary T fj also covers parts 
of interior boundary, that Helmholtz decomposition does no longer hold: There 
exist harmonic functions which are constant on different parts of T d and, hence, 
are neither in S/H], (Tl), nor in Curl Hf(Tl). 

Remark 3.6 (computation of ip). The computation of <p appears as a practical 
difficulty because <p needs to be defined through an integration of f. If f has some 
simple structure, e.g., f is polynomial, this can be done manually, while for more 
complicated f, a numerical integration of f has to be employed, but is possible in 
parallel. 

3.2 Discretization 

Let T be a regular triangulation of Tl and k G N U {0} and define 
X h (T) :=P fc (T;R 2 ) and Y h fT) := P fc+1 (T) n Y. 
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The discretization of (13.31) seeks ph. £ Xh(7) and ah £ Yh( T) with 

(Ph,qh)L 2 (Q) + (■ qh,Cuv\a h ) L 2 (n) = (p,qh)L 2 (ti) for all q h £ X ft (T), (J331a) 

(Ph, Curl P h ) L 2 (fi) = 0 for all ff h £ Y h (7). flSlJb) 

Remark 3.7. Since there are no continuity conditions on qh £ Xh(7) and since 
CurlY^(T) C Xh(7), the first equation is fulfilled in a strong form, i.e., 

p h + Curias = II k ip. 

In contrast to classical finite element methods, the approximation ph of Xu is a 
gradient only in a discrete orthogonal sense, namely (13.41 bD . For k = 0, Subsec- 
tion 1 3. 31 below proves that this discrete orthogonal gradient property is equivalent 
to being a non-conforming gradient of a Crouzeix-Raviart finite element function. 
The main motivation of the novel formulation is the generalization of this scheme 
to any polynomial degree k. 

Remark 3.8 (existence of discrete solutions). Since CurlY/,,(T) C Xh(7), the 
discrete inf-sup condition 

l|Curl/? ft || i2 ( n ) < sup ~———\\ — 1 | ^ L for all/3 h e Y h (7) 

q h &x h {7)\{0} II qh\\mn) 

is fulfilled. This and Brezzi’s splitting lemma m imply the unique existence of a 
solution ( PhiOih ) £ Xh(7) x Yh( T) to (13.41) . The equality in 

lb^llz, 2 (Q) + 11 Curl ah 11 £ 2 (Q } = ||n fe ^||^ 2{n) < |M|£a (n) . 

follows from the L 2 orthogonality of ph and Curia/!. 

The conformity of the method and the inf-sup conditions from Remarks 13.21 
and 13.81 imply the following best-approximation result. 

Theorem 3.9 (best-approximation). The solution (p,a) £ X x Y to (I3.3[) and 
the discrete solution (p/i,a/j) € X/,(7) X Yh(7) of (13.41) satisfy 

\\P-Phh^n) + ||Curl(a-a/ l )|| L2(n) (3.5) 

< ( min ||p - q h || L 2 (n) + min ||Curl(a - fi h ) || i2(n) ) • □ 

x h (T) dh& h {7) v X 

Remark 3.10. A direct analysis of the bilinear form B : (X x Yj x (I x h) -> 1 
defined by 

T,((p,a),{q,fi)) := (p,q) L 2 (f2) + {q, Curla) L 2 (n) + (p, Curl/3) L 2 (n) (3.6) 

for all p,q £ X and all a, ft € Y reveals that the inf-sup constant of B equals 5 
and, hence, the constant hidden in < in <E3 is 5. 
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Remark 3.11. The best-approximation of Theorem \3.9\ contains the term 

on the right-hand side, which depends on the choice of p. This seems to be worse 
than the best-approximation results for standard FEMs, which do not involve such 
a term. However, if p is chosen sm.ooth enough, then Curia = p — Xu has at least 
the same regularity as Xu, and therefore the convergence rate is not diminished. 
On the other hand, the approximation space for p does not have any continuity 
restriction and so the first approximation term 

min■ lb ~Qhhfin) (3-7) 

qh&x h {7) 

is superior to the best-approximation of a standard FEM, where p = Xu is approx¬ 
imated with gradients of finite element functions. However, Theorem 3.2] and 
the comparison results ofm prove equivalence of (EZD and the best-approximation 
with gradients of a standard FEM up to some multiplicative constant. 

The following lemma proves a projection property. This means that for any v G 
Hq(Q), the best-approximation of Xv in Xh( T) is a discrete orthogonal gradient in 
the sense that it is orthogonal to Curl Yj^T) and so belongs to the set of discrete 
orthogonal gradients Wh( T) defined by 

W h ( T) := {q h G X h (7) \ (q h , Curl fa) L 2 {n) = 0 for all fa G T h (T)}. (3.8) 

The projection property is the key ingredient in the optimality analysis of Sec¬ 
tion [6j 

Lemma 3.12 (projection property). It holds that Wg(7) = Hx h (T\X . More¬ 

over, if T* is an admissible refinement of T, then Rj^p-jW^T*) = Wf l (‘T). 

Proof. Let q G XHq (0). Since Curl Yh(T) C XhfiJ) and Yh( T) C Y, the orthogo¬ 
nality in the Helmholtz decomposition (13.211 implies for any fa G Yj,(T) that 

( n x h (T)9> Cur l/^)z, 2 (ft) = (q, Curl^) L 2 (n) = 0. 

This proves V//,J (Q) C Wh(7). For the converse direction, let ph G Wj l (T) 

and let u € be a solution (possibly not unique) to 

(nx h (T)^iHx h (T)Vu)L2(n) = (Phi nx h (T)Vu) L 2 (Q) for all v G Hq(Q). 

The orthogonality of pg — Hx h (T) V'u to Vi7g(H) implies the existence of a € Y 
such that ph — H^^Vu = Curia. Therefore, Curia € Xh(T) and, hence, a 
is a piecewise polynomial of degree < k + 1 and therefore a G Y/j(T). But since 
Ph G Wh(T), it holds that 

II Curl a || £ 2 (s1) = (p h ~ nx t (r)V«, Curl a) L 2 ^ = 0 

and, hence, a = 0. This proves n^mVu = Ph and, therefore, IVh (T) C 

II. v , i( , rj V//'(H). 

A similar proof applies in the discrete case and proves Hx h (y)Wh{7* = W/j(T). 

□ 
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Remark 3.13 (computational costs). Problem (13.41) is equivalent to the problem: 
Find (ph,cth) E X/^T) x Yg( T) such that 

(Curl (3 hl Curl a h ) L 2 ^ = (</?, Curl/3ft ) i2(n ) for all /3 h E Y h (T), 

Ph = - Curla^. 

Therefore, the system m.atrix is (in 2D) the same than that of a standard FEM 
(up to degrees of freedom on the boundary). 


3.3 Equivalence with Crouzeix-Raviart FEM 

The non-conforming Crouzeix-Raviart finite element space [21] reads 

ucr is continuous at midpoints of interior edges 
and vanishes at midpoints of boundary edges 


CRi(T) := | 


v cr e Pi (3") 


Since CRg(T) <2 Rq(^) the triangulation consists of more than one triangle), 
the weak gradient of a function uqr £ CRq(T) does not exist in general. However, 
the piecewise version V nc ucr E Po(T;M 2 ) defined by (V ng ucr)|t := V(wcr|t) for 
all T E T exists. The Pi non-conforming discretization of the Poisson problem 
seeks ucr E CRq(T) with 

(V nc mcr, V nc ucr)l2(q) = (/, ucr)l 2 (T 2 ) for all ucr € CRq(T). (3.9) 
The lowest-order space of Raviart-Thomas finite element functions [34] reads 


RT 0 (T) := | 


qilT E P(div, n) 


VTGT3a T GM 2 ,6 T GM 1 
with <7 rt(^) = ar + b^x ) 


(3.10) 


The Raviart-Thomas functions have the property that the integration by parts 
formula holds for functions in H( (14) as well as for functions in CRg(T). 

The following proposition proves the equivalence of the Pi non-conforming dis¬ 
cretization and the discretization (13.4p for k = 0. Note that the discretization (13.91) 
is a non-conforming discretization, while the discretization (13.41) is a conforming 
one. 


Proposition 3.14 (equivalence with CR-NCFEM). Let f G Pq(T) be piecewise 
constant and let </?rt G RTo(T) satisfy — div</jRT = f. Then the discrete solution 
( Phi a h ) E Po(T;M 2 ) X (Pi(T) D Y) to (13.4[) for k = 0 and the gradient of the 
discrete solution ucr G CRj(T) to (13.91) coincide, 

Ph = Vnc^cr- (3-11) 

Proof. The crucial point is the discrete Helmholtz decomposition [1] 

P 0 (T;M 2 ) = V NG CRj(T) © Curl(Pi(T) n Y). (3.12) 

Since p ft is L 2 orthogonal to Curl(Pi(T) fl Y), this implies ph = V nc ucr for some 
ttCR E CRq(T). Let q h = V nc ucr for some vcr G CRq(T). Then q h is L 2 
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orthogonal to Curl(Pi (T) fl Y ) and a piecewise integration by parts and (13.4|) 
imply 


(V NC ucR,V NC ncR) L 2 (n ) = (ph,Qh)L 2 (n) = (trt, qh}L 2 (n) 

= (-divt/PRT^CR^^) = (/, ^Cr)l2(q). 


Hence, ttcR = uqr solves (137911 . □ 

The projection property from Lemma 13.121 generalizes the famous integral mean 
property Vnc-ZncT = npo( 0 ';R 2 )^ 7 ' u for all v € Hq(FI) of the non-conforming interpo¬ 
lation operator 1^. 

Remark 3.15 (higher polynomial degrees). For higher polynomial degrees k > 1. 
the discretization (E3D is not equivalent to known non-conforming schemes \2B, 
[20 i '. 21\ 1 28fi in the sense that Wh( T) V NC V/ l (T) for those non-conforming finite 
element spaces V/ l (T). This follows from V NC V/ l (T) ^ W/j(T) for non-conforming 
FEMs with enrichment. A dimension argument shows 

dimmer)) > dim 14 (“J) 

for the non-conforming FEMs of [Uffl/ without enrichment and therefore Wh( T) 7 ^ 
V nc 14(T). Moreover, this proves that the generalization of the projection property 
to higher m.oments from. Lemma 1 3.12\ cannot hold for those finite element spaces, 
in contrast to the discretization (ED). 

4 Extensions 

Subsection 14.11 generalizes the novel FEM to quadrilateral meshes and proves a new 
discrete Helmholtz decomposition for the Q 1 rotated non-conforming Rannacher- 
Turek FEM |33 j. Subsection 14.21 discusses a discretization with Raviart-Thomas 
functions. 

4.1 Quadrilateral finite elements 

For this subsection, consider a regular partition 7 of H in quadrilaterals. Define 
for the reference rectangle T = [0, l] 2 

Qk(T) ■■= {v h <E P- 2 k{T) | 3/,p <E Pk([0, 1]) : v h (x,y) = f(x)g(y)}. 
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Given TgT, let ifx ■ T —>• T denote the bilinear transformation from the reference 
rectangle to T. For consistency, let P_i([0,1]) := {0} and set 


V Q , k ( 7 ) := {/3j 

xr x (T) ■■= 


VT € 7 : (Ph°i’T)\f£Qk(T)\, 
3a E R, 6, c 

r h EL 2 (T;R 2 ) 


== 


T t a 2 (!l;R 2 ) 


3a E R, b,c E P fe _ 2 ([0,1]), d, e E Q k -i(T) 
such that V(x,y) E T 

fx k b(y) + d(x,y) 
\y k c(x) + e(x, y) 

VP E 73pT E X[ ect (f) such that 

(pi 0 Vr)lf 


Th(x,y) = a 


-x k y k 1N 
x k ~i y k ~*~ 


0 1 

-1 0 


Dty, 


— 1\T 


1 IpT 


0 -1 

1 0 


Pt 


Then a discretization with respect to the quadrilateral partition seeks p k E _A£ ect (T) 
and a h E Vq^(T) with 

{ph,qh)LZ(n) + {qhi Curl«/ i)l 2 ( q) = (<p,q h )L2(n) for all qh e ^T ct (T), 

(Ph, Curl/ 3 / l ) L 2 ( Q) = 0 for all f3 h E Vq^T). 

Let /3/j E Vq^T), i.e., (/H^ ° V’lOI'f E Qk(T). A direct calculation reveals for all 
T E T 

' 0 1 
-1 0 

0 1 
-1 0 


((Curl /3 h ) T oip T )\ T = (V^oi/jyo ^ 1 )) 1 oip T 
= (X(/3 h o' l f T )) T Dii/j- 1 ) o-0r 


Let (/^o i/j T )(x,y) = (ax k + f{x)){by k + g(y)) with a, b E R and /,5 E P fc _i([0,1]). 
Then it holds 


V(/3/ l o if T ) = ab/c I Zfc-i^k + 


x k 1 jf k 


y x 


by k df(x)/dx\ f akx k 1 g{y) + g{y)d f (x) / dx\ 


+ 


a x k dg(y)/dyj \bky k f(x) + f(x)dg(y)/&y J 


and therefore 
(V(/3 h o^t)) T 


0 -1 

1 0 


= nbfc 

V V-: 

=: ( p T (x,y)) T 


y k 1 x k 
-x k ~ 1 jf k ~*~ 


ax k dg{y)/dy 
— by k df{x)/dx^ 


bky k x f{x) + f(x)dg{y)/dy \ 
-akx^gify) - g(y)df (x)/dxj 


This implies pr E Xjj ect (T) and (V(/3/ l o V’t)) = (0,1; —1,0 )pr- The combination 
of the previous equalities leads to 


((Curl j3 h ) o V>t)|t = 


_°1 \) D ^t 1 ) T °^T 


0 -1 

1 0 


Pt- 


Consequently, Curl E A[ cct (T). This and the conformity of the method prove 
as in Section [3] the following statements 
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(i) unique existence of solutions, 

(ii) the best-approximation result 

\\P ~ Ph\\tf(ci) + ||Curl(a-a ft )|| L 2 (n) 

~ l|p “ qhU ^ + fctm l|Curl( “ ' A)lli> «“>)' 


(iii) the projection property 

n^ect^vtfo 1 ^) c wr\7) 

for 

wr\7) = {q h G X£ ect (T) I V& G V Qtk (7) : (q h , Curl&) La(n) = 0}. (4.1) 

Remark 4.1. The properties (i)-(iii) still hold for any XhifJ) with X(. ect (7) C 
X h {7) C X. 


The remaining part of this subsection proves the equivalence of the lowest-order 
rectangular discretization with the non-conforming Rannacher-Turek FEM |33| . 
To this end, define for the reference rectangle T and the bilinear transformation 
-0T : T -> T, 


Q m \T) : = 

^ 0t (T) : = 


span{l ,x,y,x 2 - y 2 }, 


v h G L 2 (Ll) 


VT G 7: (v h o if T )\f G Q rot (f) and 
f E Vh ds is continuous for all interior 
edges E and vanishes at boundary edges E 


(4.2) 


The following lemma proves a relation between the cardinalities of the quadri¬ 
laterals, nodes, and interior edges of a quadrilateral partition similar to Euler’s 
formulae 

card(S) + card(£(f!)) = 3card(T), 
card(£(0)) + card(N) = 2card(T) + 1 

on triangles. This enables a dimension argument in the proof of the discrete 
Helmholtz decomposition in Theorem 14.31 below. 

Lemma 4.2 (Euler formula for quadrilateral partitions). Let 7 be a regular parti¬ 
tion of LI in quadrilaterals with edges 8, interior edges 8(LI), and vertices N. Then 
it holds that 3card(T) + 1 = card(£(fl)) + card(N). 

Proof. Define a triangulation Ta of fl in triangles by the division of each quadri¬ 
lateral into two triangles by a diagonal cut. Let £a denote the edges of Ta, £a(^) 
the interior edges and Na the vertices. Then the following relations between the 
two partitions hold 


card(TA) = 2card(T), card(£A) = card(£) + card(T), 

card(£A(f^)) = card(£(D)) + card(T), card(NA) = card(iH). 
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This and Euler’s formulae for triangles (14.31) prove 

card(£(fi)) + card(AT) = card(£A(^)) — card(T) + card(NA) 

= 2card(‘J'A) + 1 — card(T) = 3card(T) + 1. □ 

The following theorem proves that the solution space Wfi ct (7) from (14.ip 
equals the piecewise gradients of functions in on a partition in squares 

for k = 1. 


Theorem 4.3 (discrete Helmholtz decomposition on squares). Let 7 be a regular 
partition of fl in squares. Then, 


Xl ec \7) = V NC 4^ ot (T) © Curl Vgq(T) 


and the decomposition is L 2 orthogonal. 


(4.4) 


Remark 4.4. The L 2 -orthogonality in (ED still holds for a partition in parallel¬ 
ograms. However, V NC LJJ ot (T) % X[ ect (T) for general quadrilateral partitions. 


Proof of Theorem \4.S\ Let Vh £ L^c^T) and fih £ Vq i i(‘T). A piecewise integra¬ 
tion by parts leads to 


(Vnc!^, Curl^)^2(Q) 


/ [vhlEVPh ■ t e ds. 
egE J E 


Since 7 consists of parallelograms, the bilinear transformation fix ■ T —> T is 
affine and, hence, fih\E is affine on each edge E £ £. This implies that V/3/j • 
te is constant. Since the integral mean of \v^\e vanishes, this proves the L 2 
orthogonality. 

Let Vh £ L(™ t (T). A computation reveals for all T £ 7 that there exist /r£l 
and gx £ M 2 such that 


Vv h (x,y) = D(fi T 1 ) T 





For k = 1, X{ ect (T) reads 




< T h £ L 2 (H;M 2 ) 


VT £ 73ax £M.,dx£ M 2 such that 

(■th ° V’t) if = ^ ^ 



Since all T £ 7 are squares, Dfix and ( 0 , 1 ;— 1 , 0 ) commute, and, hence, Vu/, £ 
Xf ct (T). Thus, V N c^ ot (T) 0 CurlF Q) i(T) C Xf ect (T). The dimension of 
equals card(£(H)) and the dimension of Curl Vq i(T) equals card(Af) — 1, while the 
dimension of Xf ect (T) equals 3card(T). This and Lemma 14.21 prove the asser¬ 
tion. □ 


Remark 4.5 (arbitrary quadrilaterals). The best-approximation (ii) from above 
proves quasi-optimal convergence even for arbitrary quadrilaterals. Standard in¬ 
terpolation error estimates for Vq i(T) and for Pq{7'M. 2 ) C X[ ect (T) \19j lead to 
first-order convergence rates of h for sufficiently smooth solutions. This should 
be contrasted with where quasi-optimal convergence is only obtained for a 

modification of (14.21) where Lj^ )t (T) is defined in terms of local coordinates. 
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4.2 Relation to mixed Raviart-Thomas FEM 


This subsection shows that the classical mixed Raviart-Thomas FEM [34| can be 
regarded as a particular choice of the ansatz spaces in the new mixed scheme. 

Let 7 denote a regular triangulation of hi in triangles. Define the space of 
Raviart-Thomas functions [34] 

*Ht(T) = Urt G H (div, fi) I VT E 7 : q RT \ T (x ) € R fc (T; M 2 ) + P fe (T) x} 

and 


T rt (T) := P k+1 (7) n Y. 


Then the following problem is a discretization of (13.31) : Seek (prt, ort) € A’rt(T) x 
Yrt(T) with 


(prt,9rt)l 2 (Q) + (<7RT, CurlaRT)L2 ( n) = (<P,Qrt) 
(prt , Curl /3 rt ) tP (O) = 0 


for all gRT E 3 lrt(T), 

(4.5) 

for all /3 rt G Trt(T). 


Since Curl Trt(T) C Pfc(T;M 2 ) and div Curl urt = 0 for all urt € Yrt(T), it 
follows Curl Irt (T) C Art (T). This and the conformity of the method guarantee 
as in Section [3] and in Subsection 14.11 the unique existence of solutions, a best- 
approximation result, and the projection property 


n.Y RT (T)ViLo(D) C VFrt(T) 

:= {<?RT G A'rt(T) | V/3rt G Trt(T) : (^rt, Cuii/3Rr)z, 2 (n) = 0}. 

The discrete Helmholtz decomposition of mum proves 

Xrt(T) = V RT P k (7) ® Curiy RT (T) 

with the operator Vrt : P k {7) Art(T) dehned for all urt G P k (7) by 

(Vrturt, (Zrt)l 2 (Q) = — (^rt j div g R x ) l 2 (Q) f° r all Qrt G Ar T (T). 

This decomposition yields the equivalence of (14.51) with the problem: Seek (prt 5 urt) G 
Ar T (T) x P k (7) with 


Prt = Vrturt, 

(wrt, divp RT ) L 2 (n) = (div n XRT(T) <£, wrt)l2 (H) for all ^RT G P k {7). 

This is the classical Raviart-Thomas discretization with / replaced by div nX rtPRP- 
Assume now that the right-hand side ip G Art(T) is a Raviart-Thomas func¬ 
tion. Since by definition Yrt(T) = Yjj(T) with Yh(7) from Subsection 13.21 and 
since «rt is the solution of 


(Curl/3 R T, CurlaR T ) L 2 (f 2 ) = (<p, Curl/3R T ) L 2 (n) for all /3 RT G Fr T (T), 

it holds ckrt = ah with ah from (13.41) . Since ip = prt + CuHort and H Xh pj^ip = 
Ph + Curias, it follows 


Ph = n^(T)PRT- 

For k = 0, the equivalence with the Crouzeix-Raviart FEM (13.111) then proves the 
identity 


Vng'Ucr = n Xh ( T )PRT, 

which is also known as Marini identity 13 EZj. 
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5 Medius analysis 

The medius analysis of 1251 H5] proves for the discrete solution ucr CRq(T) to 
(13.91) the best-approximation result 

||Vnc(u - u C R,)||L2(n) ^ min ||V NC (u - u C r)||l 2 (0) + osc(/, T). (5.1) 

^crgCRI(O-) 

The following theorem proves a generalization for the discretization (13.41) for the 
lowest order case k = 0 . 

Theorem 5.1 (best-approximation property). Let ( p,a ) £ X x Y be the solution 
to (13.31) and (ph,cxh) € Tb(T;IR 2 ) x (Pi(T) fib) be the solution to (13.4p . Then the 
following best-approximation result holds 

I \P ~ Phh 2 (n) ^ lb - noplUyo) + osc(/, T) 

gu (/,^Cr)l2(Q) ~ (<P, VNc^CR)L2(q) (5.2) 

t >CR€CRj(T)\{0} ||V NC UCR||L2(n) 

Remark 5.2. If ip is a lowest-order Raviart-Thomas function, then it allows for 
an integration by parts formula also with Crouzeix-Raviart functions (see Subsec- 
tion \3.3\) . Therefore, the third term on the right-hand side of (I5.2[) vanishes. This 
and the equivalence with the non-conforming FEM of Crouzeix and Raviart from 
Subsection 1 3 . 31 reveal the best-approximation result m- 

The remaining part of this section is devoted to the proof of Theorem 15.11 
The following lemma from na ei is the key ingredient of this proof. Recall the 
definition of CRj(T) from Subsection 13.31 

Lemma 5.3 (companion). For any ucr £ CRq(T) there exists v £ Hq(£ 1) with 
the following properties 

(i) n 0 V NC (n - v C r) = 0 , 

(ii) n 0 (u - u C r) = o, 

(Hi) Ibj^cR - v)\\ L 2(n) + ||V N c(ucr - v)\\ L 2(ty < HVncUcrIIl 2 ^)- □ 

Proof of Theorem \5. 1\ Define qh := Uop—ph £ Po(T; R 2 )- The projection property 
of Lemma [3T2] implies that q/ 1 £ Wh(7) and the discrete Helmholtz decomposition 
(13.121) guarantees the existence of ucr £ CRg(T) with qh = V nc ucr- 
Hq(Q) denote the companion of ucr from Lemma [5.31 Then 

(p ~ PhiQh)L 2 (n) = (p, V nc (u C r - v)) L 2 ( n) + (p, Vu) L 2 (0) 

- (p/i, Vnc^Cr)l 2 ( 0 )- 

The properties (i) and (iii) from Lemma 15.31 yield for the first term on the right- 
hand side 


Let v £ 

(5.3) 


ip, Vnc(^cr - (i) (ii) * * v ))l 2 (ci) = (p- n 0 p, V nc (u C r - v)) L 2 (n) 

^ lb - n oP||L2(n) ||V NC ucR||L2 ( n). 


(5.4) 
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The problems (13.3D and (13.4D lead for the second and third term on the right-hand 
side of (15.3D to 

(p, Vv) i 2(Q) - (p h , V NC ncR)L2(Q) = (p, Vn) L 2(Q) - (ip, V NG «CR)L 2 (n)- 
Since — div<£> = /, it follows 

(<P,Vv) L 2(f2) - (<p, V NC fcR)i,2(n) 

= (f,v -vcr)l 2 (q.) + (/^cr)l 2 (q) - (<£, V nc ?;cr)l 2 (q). 
Properties (ii) and (iii) of Lemma 15.31 prove 

= (f,v-v cr)l 2 (q) <osc(/,T)||V NC u C R||L 2 (n)- 

The combination with (15.3|> and (15.41) and a Cauchy inequality yield 


(P ~ Ph,qh) L 2 (H) ^ \ \\P ~ n 0 p||L 2 (n) + osc(/, T) 

(f,vcR)L 2 (n) ~ (<P, V nc u C r)l2(q) 


+ sup 
^c R eCRj(T)\{0} 


II V N c^cr||z,2( 0 ) 


11% Hl 2 (Q) ■ 


This and 

b-Pfcllia(n) = lb- n oPllL2(o) + II%IIl2(o) = lb- n oPll|2 (n) + (p-Ph,qh)iP{a) 

prove the assertion. □ 

Remark 5.4 (higher polynomial degrees). For k > 1, Remark I ff. 4 51 implies that 
an analogue of Lemma 1 5. .4 cannot be proved in the same way. 

6 Adaptive algorithm 

This section defines an adaptive algorithm based on separate marking and proves 
its quasi-optimal convergence. 

6.1 Adaptive algorithm and optimal convergence rates 

Let To denote some initial shape-regular triangulation of 0, such that each triangle 
T £ T is equipped with a refinement edge Et £ £(T). A proper choice of these 
refinement edges guarantees an overhead control [7]. 

Let T (N) denote the subset of T of all admissible triangulations with at most 
card(To) + N triangles. The adaptive algorithm involves the overlay of two admis¬ 
sible triangulations T, T* £ T. which reads 

T <gi T* := {T G T U T* | 3K G T, K\ G T* with T CKF A*}. (6.1) 
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Given a triangulation T#, define for all T £ the local error estimator contri¬ 
butions by 


A 2 (T,,T) := \\h 7 cui] NC p h \\ 2 L2{T) + hT Y W\Ph\E ■ te\\ 2 L 2 (e) , 

egE(t) ( 6 . 2 ) 



p 2 (T) 

:= \\<p 

^-kT Hl 2 (T) 


and the global error 

estimators by 


A? 

:= A 2 (T,,T,) 

with 

A 2 (T £ ,M) := 

Y A 2 (T £ , T) for any M C T e 





TeM 

Fi 

:= 

with 

/r 2 (M) := 

Y^ h 2 {T) for any MCTj 





TeM 


The adaptive algorithm is driven by these two error estimators and runs the fol¬ 
lowing loop. 

Algorithm 6.1 (AFEM). 

Input: Initial triangulation To, parameters 0 < 6 a < 1, 0 < ps < 1, 0 < k. 

for i = 0,1, 2, ... do 

Solve. Compute solution (pi,o.e) £ Xh(Ti) x Y/ l ( 7e) of 113.411 with respect 
to T£. 

Estimate. Compute local contributions of the error estimators (A 2 (T^, T)) Tg ^ 
and (ii 2 {T)) Te7e . 
if fij < nX 2 then 

Mark. The Dorfler marking chooses a minimal subset JvQ C Tf such that 

e A \j < a 

Refine. Generate the smallest admissible refinement T^ + i of T^ in which 
at least all triangles in are refined. 

else 

Mark. Compute a triangulation TgT with /i 2 (T) < p_B/U 2 . 

Refine. Generate the overlay T^+i of and T. 

end if 
end for 

Output: Sequence of triangulations (T^ eNo , discrete solutions (p^a^eNo and 
error estimators (A^)^ g pj 0 and ♦ 

Remark 6.2 (separate versus collective marking). The residual-based, error esti¬ 
mator ^ X 2 + p 2 involves the term ||y> — without a multiplicative positive 

power of the mesh-size. Therefore, the optimality of an adaptive algorithm, based on 
collective marking (that is k = oo and A replaced by \J A 2 + p 2 in Alaorithm \6.1\ ) 
does not follow from, the abstract framework from m The reduction property 
(axiom (A2) from JT3/), is not fulfilled. Alaorit,hm \6.1\ considered here is based on 
separate marking. In this context, the optimality of the adaptive algorithm (see 
Theorem m can be proved with a reduction property that only considers A. 

Remark 6.3. The step Mark in the second case (p 2 > kX 2 ) can, be realized by the 
algorithm Approx from Um, i-e., the thresholding second algorithm followed 
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by a completion algorithm. For this algorithm, the assumption (Bl) optimal data 
approximation, which is assumed to hold in the following, follows from the axioms 
(B2) and (SA) from Subsection \6.5\ Jffl j . For a discussion about other algorithms 
that realize Mark in the second case, see m- 

For s > 0 and (p,a,ip) £ X x Y x H(div, fi) define 
I(P>“>¥0L. : = SU P 

NeNo JfcJipV) \ 

+ Jnf l|Curl(a - A-)|| La(n) + \\ip - n Xfc(T) ¥?|| £ 2 (n) ) . 

Remark 6.4 (pure local approximation class). Since 0 is assumed to be a Lip- 
schitz domain, all patches in an admissible triangulation TgT are edge-connected, 
i.e., for all vertices z £ IN’ and triangles T,K £ T with z £ T (~l K, there exists 
m £ No and Kq, ..., K m £ T with Kq = T, K m = K, z £ Kq D ■ ■ ■ D I\ m and 
Kj—iDKj £ £ for all 1 < j < m. Under this assumption, \f2\ Theorem 3.2] shows 

min \\V(v - v h )\\ L 2 ^ & ||Vu - n fe Vv|| £ 2 (n) for all v £ H 1 (Q). 
v h GP k+1 (7)nHUn) v ; v ; 

Hence, 

\{P^,g>)\ As « |(p,a,<p)U' 

:= S ) JV ‘d}( f i v)( l|j, “ n - v ''< T > p|lLain > 

+ 11 Curl a - n x ^ (T) Curla|| L2(Q) + ||p - n Xh(T) ^|| £ 2 (n) ) . 

In the following, we assume that the following assumption (Bl) holds for the 
algorithm used in the step Mark for fi] > kX] (see Remark 16.31) . 

Assumption 6.5 ((Bl) optimal data approximation). Assume that 

is finite. Given a tolerance Tol, the algorithm used in Mark in the second case 

(/ij > K.Xj) in Alaorithm \6.1\ computes T* £ T with 

card(T*) — card(To) < To 1 -1 ^ 2ct ) and /i 2 (T*) < Tol. 

The following theorem states optimal convergence rates of Algorithm 16. II 

Theorem 6.6 (optimal convergence rates of AFEM). For 0 < pb < 1 and suf¬ 
ficiently small 0 < k and 0 < 0 < 1, Algorithm 16.11 computes sequences of tri¬ 
angulations (T^gN and discrete solutions (pt, ag )^gN f or the right-hand side ip of 
optimal rate of convergence in the sense that 

(card(T^) - card(T 0 )) a (||p - p*|| £ 2 (n) + ||Curl(a - a^)|| £2(n) ) < \(p,a,ip)\ As . 

The proof follows from the abstract framework of [16], which employs the 
bounded overhead [7] of the newest-vertex bisection, under the assumptions (Al)- 
(A4) and (B2) and (SA) which are proved in Subsections 16.2116.51 
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6.2 (Al) stability and (A2) reduction 

The following two theorems follow from the structure of A. 

Theorem 6.7 (stability). Let T* be an admissible refinement of 7 and M C TflT*. 
Let (py*,ay*) S A ft (T*) x 1^(7*) and {py,ay) £ X h (7 ) x Y h {7) be the respective 
discrete solutions to (13.41) . Then, 

|A(T*, M) - A(7,M)| < |K -p T || L 2 (n) . 

Proof. This follows with triangle inequalities, inverse inequalities and the trace 
inequality from [TOl . p. 282] as in [181 Proposition 3.3]. □ 

Theorem 6.8 (reduction). Let 7* be an admissible refinement of 7. Then there 
exists 0 < p 2 < 1 and A 2 < 00 such that 

A 2 (T*, 7* \ T) < p 2 A 2 (T,T\T*) + A 2 |K -Pr|li a(n) - 

Proof. This follows with a triangle inequality and the mesh-size reduction property 
/i.jjy < |r/2 for all T £ 7* \ 7 as in [18, Corollary 3.4]. □ 

6.3 (A4) discrete reliability 

The following theorem proves discrete reliability, i.e., the difference between two 
discrete solutions is bounded by the error estimators on refined triangles only. 

Theorem 6.9 (discrete reliability). Let 7+ be an admissible refinement of 7 with 
respective discrete solutions (py*,oty*) £ A^(T*) x 4^(7*) and (py, ay) £ A'^(T) X 
Y h (7). Then, 

IIpt -PrJhm + IICurlK - ao-Jll^n) < A 2 (7,7\7*) + /r 2 (7,7\ 7*). 

Proof. Recall the dehnition of Wh(7 .*.) from (13.81) . Since py — py* £ X/,(T*), there 
exist ay* £ W h (7fi) and ry* € Y h (7±) with py - py* = ay* + Curl ry*. Since 
W h { T*)T L2(n) Curiy fc (T*), 

lkT*lli 2 ( 0 ) + II Cui ' 17 XIIl 2 (Q) = l|PT-PT*|li 2 (n) - 

The orthogonality furthermore implies that the discrete error can be split as 

l|PT-PT*||£ 2 (n) = (PT-PT*,0- T j £ a(n) + (py - Py*,Curlry*) L 2 {n) . 

The projection property, Lemma T3.121 proves 11^mo'er* £ 1T/j(7). Hence, prob¬ 
lem (13.41) implies that the hrst term of the right-hand side equals 

(pr -PT*,crT*)L 2 (n) = {^X h (7)T ~ </?> )l 2 (H) = {^-X h (7)P - KlL) 1 / 3 ! ■ 

For any triangle T £7 D 7*, it holds (n^m^ — Y Xh (y*)< p)\t = 0. Therefore, 

(Rx h {T)V ~ n Yh ( 0 '*) < / :, ’ cr a'jL 2 (a) < l|nx h (7)¥ , -nx h (T*) ( PllT\T* IKIKn). 
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Since T* is a refinement of T, it holds 

linx^T)^ -- p)IIt\t* < llp-nx^pIlTyr*- 

Let r 7 G Yh(7) denote the quasi interpolant from [35] of r 7t which satisfies the 
approximation and stability properties 

II^Vt* - ri)\\ L 2 {Q) + ||Curl(r T * - r 7 ) || £2(n) < ||Curlr T *|| L 2 (n) 

and (rj) | e = (^t*)|_b for all edges E G £(T) n £(T*). Since pj G Wh(7) and 
Pt* G W h ( T*), 

(pt - Pt*, Curlr T J i2(r2) = (p T , Curl (r 7t - r T )) L 2 (n) . (6.4) 

An integration by parts leads to 

(PT, Curl(r T * - r T )) i2(n) = -(curl NC p T , r 7ir - r 7 ) L 2 (a) 

+ ^2 / br ' T £k( r T* ~ Cr) ds - 

Te£(T) Je 

For a triangle T G T fl T*, any edge E G £(T) satisfies A G £(T) n £(T*). Hence, 
( r T)|r = ( r T*) |t f° r all T G 7 nT*. This, the Cauchy inequality and the approxi¬ 
mation and stability properties of the quasi interpolant lead to 

-(cuHnc^tx - r 7 ) L 2 ( 0 ) % ||/lo-cur^c po-ll^ ||Curlr T J| L2(n) . 

Since (r 7 )\ E = ('/"X)| E for all edges E G £(T) fl £(“X*), the approximation and 
stability properties of the quasi interpolant and the trace inequality [10, p. 282] 
lead to 

V / [pt • r E \ E {r 7ir - r 7 ) ds 

e^ Je , 

I -^- ( 6 - 5 ) 

SJ 22 h T \\\p7 ■ te\e\\ 2 L 2 {e) ||Curl7x|| £ 2 (n) - 

Y £e£(T)\£(X) 

The combination of the previous displayed inequalities yields 

IIpt - Pt* ||£ a(n) < A 2 (T, 7 \ T*) + /r 2 (T, T \ T*). 

Since Curl a 7 = X\ Xh ( 7 )T - Pt and Curias* = n Xfe(T ^ - jx*, the triangle in¬ 
equality yields the assertion. □ 

The discrete reliability of Theorem 16.91 together with the convergence of the 
discretization proves reliability of the residual-based error estimator. This is sum¬ 
marized in the following proposition. 

Proposition 6.10 (efficiency and reliability of the residual-based error estimator). 
Let (p,a) G X x Y and (ph,cth) G Xh(7) x Y),(T) be the solutions to (13.31) and 
(El for some TgT. There exist constants C e s,C Te \ > 0 with 

+ P 2 (^)) - IIP-PfellL 2 (n) + l|Curl(a - a/Olli^n) 
<C 2 el (A 2 (T,T)+/r 2 (T)). 

Proof. The a priori error estimate from Theorem 13.91 implies the convergence of 
the discrete solutions. This and Theorem 16.91 proves the reliability. The efficiency 
follows from the standard bubble function technique [43] . □ 
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6.4 (A3) quasi-orthogonality 

The following theorem proves quasi-orthogonality of the discretization (13.4|) . 

Theorem 6.11 (general quasi-orthogonality). Let (7j \ j £ N) be some sequence 
of triangulations with discrete solutions ( pj,aj) £ Xh (T j ) x Yg(flj ) to (13.411 . Let 
t £ N. Then, 

OO 

y: (ll Pj ~Pj- iI!l 2 (o) + ||Curl(c>!j - Oij-i)\\ 2 L 2 (ty ) < A|_! + /x|_i. 
j=t 

Proof. The projection property, Lemma 1.3. 1 21 proves G Wh(fTj-i) with 

Wh(7j-i) from (13.811 . Hence, problem (13.411 leads to 

{Pj-i,Pj - Pj-l)L 2 (Q .) = Xhpj-^Pj -Pj-i)L 2 (ty, 

(PjiPj ~Pj-i)L 2 (n) = (P>Pj) - (Vi^-x h (7 j _ 1 )Pj') l 2 (Q)- 

The subtraction of these two equations and an index shift leads, for any M £ N 
with M > £, to 

M M 

Y ii pj -pj- iiii 2 (o) = Y(v>Pj - n 

j=e j=i 

M M -1 

-Y& ^XhiTj-YP^L 2 ^) + Ys (PiPj)L 2 (n) (6.6) 

j=e j=t —i 

M 

= (ViPt-i Pm)L 2 (fi) + 2 YXt-Pj ~ n x,4T 3 _i)Tj)L 2 (n)- 

j=e 

Since — If x h (7--i)Pj £ ^(Tj) is L 2 -orthogonal to X^Tj-i), a Cauchy and a 
weighted Young inequality imply 

M 

2 Y^ ,P ^ ~ n XhffTj^Pj)^^) 
j=l 

M 

= 2 Y^X h (7j)P ~ n XhVj-^ViPj ~ Ha - h (7 j _ 1 )Pj)L 2 <si) (6.7) 

3=t 

M j M 

^ 2 Y \\ u x h (7,)<P - n^ (T ._ 1 )^||i 2(n) + - Y II Pj ~ nxhfr-jPjWhw 

j=t j=t 

The orthogonality Ux^ip - U Xh(: j j _ rn )P±L 2 (n) x h(7j- m ) for all 0 < m < j 
proves 

M 

Y IIHatCYT’ - n A fe (T J _ 1 )7 , ll|2(n) = Iln^^M)^ - n x h (^_ 1 )¥’lli2(n)> ( 6 -8) 

j=i 
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The definition of fig yields 

ll n -Y h (T M )^ - n^cr^o^ll l 2 (Q) = ll n x h (T M )( < / ? - ^ 

< Ft-\ ■ 

The combination of (ESD (EH and ||pj - n^(T^^'II l 2 ^) < ||pj - Pj-ih^Q) 
leads to 

1 M 

2 S HPi _ ^-illl 2 (n) < 2 l4 -1 + {p>,Pi-i~Pm)l 2 {Q)- (6.10) 

j=i 

The combination of the arguments of (16.4D — (16.51) proves 

(Curl(a M - a£_i),p£_i) L 2 (n) < Af_i ||Curl(a M - &g-i)\\ L 2 (n) ( 6 . 11 ) 

This, the discrete problem (13.4| . and the discrete reliability ||Curl(aA/ — ag-\)\\ L 2 ^ 
\g-i + Hi -1 from Theorem 16.91 lead to 

(pi-i ~ PmL 2 (p.) = {Pt-i ~ Pm,P i-i + Curlar_i) L 2( 0 ) 

= ( Pe-i ~ PM,Pi-i)L 2 (fi) = (Curl(aM — ag-i),pg-i) L 2 ^ 

^ A t-i ||Curl(aM — a t-i)llL 2 (n) ~ ('V-i + Pi-i) 2 - 

This and a further application of Theorem 16.91 leads to 
{p,Pe-i ~ PM)L 2 (n) 

= {ip - n x h (7 t _ 1 ) l P,P(.-l - PM)L 2 (fl) + {pi -1 - PM,^X h (7 e _ 1 )P)L 2 (0.) 

< \\ip - n Xh (T £ _ 1 )¥ 5 ll l 2 (o) \\pe-i ~ Pm\\l 2 (q) + (A«-i + Pt-i) 2 L 2 (si) 

< {Xg-i + pg-i ) 2 ■ 

(6.12) 

The combination of (16.101) with (I6.12|) implies 

M 

Y. II Pj ~Pj- illi 2 (o) ~ ^e-i +3-1- (6.13) 

3=1 

The Young inequality, the triangle inequality, and Curl aj = Ux h ( t T ] ) ! p ~Pj imply 
M 

5^||Curl(aj -aj-i)||i2 (n) 

3=1 

M M 

- 2 Y ll p J ~ P3-t\\ 2 L 2 (p) + 2 Yj ll n -Y(T')^ " n x h (Tj_i)Vll|2(0)- 
j=i j=e 


Since M > f is arbitrary, the combination with (16.8|) . (16.91) . and (16.131) yields the 
assertion. □ 
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6.5 (B) data approximation 

The following theorem together with Assumption 16.51 form the axiom (B) from [T 6 j. 

Theorem 6.12 ((B2) quasimonotonicity and (SA) sub-additivity). Any admissi¬ 
ble refinement T* of T satisfies 

M 2 (T*) < /t 2 (T) and ^ ^ 2 ( T ) ^ I j2 ( K ) f° r al1 K G T ' 

Te T* 

TCK 

Proof. This follows directly from the definition of /a. □ 

7 Extension to 3D 

This section is devoted to the generalization to 3D. Subsection 17.1 I dehnes the novel 
discretization and comments on basic properties, while Subsection 17.21 is devoted 
to optimal convergence rates for the adaptive algorithm. 

7.1 Weak formulation and discretization 

For this section, let SI C K 3 be a simply connected, bounded, polygonal Lipschitz 
domain in R 3 . For the sake of simplicity, we also assume that <9fl is connected 
(i.e., fl is contractible). The Curl operator acts on a sufficiently smooth vector 
held /3 : SI —> R 3 as Curl/3 = V A ft with the cross product or vector product A. 
Let I7(Curl,fl) denote the space of all /3 E L 2 (fl;R 3 ) with Curl/3 E L 2 (fl;R 3 ) for 
the weak Curl, i.e., 

f v-Cur\/3dx= f /3-Curlvdx for all v E C"?°(fl; R 3 ). 

J Qj J 

In contrast to the two-dimensional case, iL(Curl, fl) ^ i7 1 (D;M 3 ). The Helmholtz 
decomposition in 3D reads 

L 2 (fl;R 3 ) = Vflo( n ) © Curl (Curl, fi) 

and the sum is L 2 orthogonal. It is a consequence of the identity 

{r E H (div, fl) | div r = 0} = Curl If (Curl, fl) 

in the De Rham complex [9]. 

Let ip E H( div, fl) with — div</? = /. Then the Poisson problem 
equivalent to the problem: Find (p,a) E L 2 (fl;R 3 ) x L7(Curl, fl) with 

(P,q)l2(Q) + (< 7 , Curla)^ 2 (Q) = (<p,q) L 2 (n) for all q E L 2 (D;R 3 ), 

(p, Curl/ 3 ) L 2 (Q) = 0 for all /3 E H(Cm\, fl). 

In contrast to the two-dimensional case, the operator Curl : H (Curl, fl) — > L 2 (fl; R 3 ) 
has a non-trivial kernel. Classical results |35| characterize this kernel as ViL 1 (fl). 


(7.1) 

(13.11) is 

(7.2) 
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To enforce uniqueness, we can reformulate m as follows. Seek (p, a , w) G 
L 2 (ft;R 3 ) x L7(Curl, 17) x (iL^ft) D Lg(ft)) with 

(P,Q)L2(ci) + (g,Curla) i2 ( n) = (ip,q) L 2 ( n) for all q G L 2 (ft;R 3 ), 

(p, Curl/3) L 2 (r2) + (/?, Vw) L 2 (q) = 0 for all P € if (Curl, 0), 

(a, V»)i 2 (n) = 0 for all v G (if 1 (ft) D ig(ft)). 

Note that {/? G if (Curl, ft) | Curl fj = 0} = Vif 3 (ft) implies w = 0. 

Standard finite element spaces to discretize H (Curl, ft) in 3D are the Nedelec 
finite element spaces [501,131] (also called edge elements) which are known from the 
context of Maxwell’s equations. Let T be a regular triangulation of ft in tetrahedra 
in the sense of [19]. The spaces of first kind Nedelec finite elements read 

Y Nik (T) := P k (T ; R 3 ) + (iA P k (T ; R 3 )), 

Y N , k {7) := {p h € if (Curl, ft) | VT G T : p h \ T G Y N>k (T)}. 

Let X k (7) := P k ( T;R 3 ). Since CurlYr^fcOO C X k (7), a generalization of (13.41) to 
3D seeks (ph,ath) G X h {7) x lN,fe(T) with 

(Ph,Qh)L*{n) + (^,Curla ft ) L 2 (n) = {<p,qh)L 2 (Sl) for all q h G X h (7), 

(. Ph, Curl Ph) L 2 (fi) = 0 for all p h G T Njfc (T). 

The discrete exact sequence [9] implies that the elements in Y^ jk (7) with vanishing 
Curl are exactly the gradients of functions in U k {7) := P k +i(7) D ff 1 (ft) D L 2 (ft). 
Therefore, the uniqueness in (17.3p can be obtained in the following formulation. 
Seek (ph, a.h-,Wh) G X h (7) x Yn, k (7) x U h (7) with 

(Ph,qh)L 2 ( n) + {qh, Curla h ) L 2 (n) = {v,Qh)L 2 (n) for all Qh G X h {7), 

{Ph, Curl/3 h ) £ ,2 ( Q) + {P h , Xw h ) L 2(n) = 0 for all p h G In./cC^), (7.4) 

(a h , Vv h ) L 2(n) = 0 for all v h G U h {7). 

Note that VC/^(T) is the kernel of Curl : ^^(T) —> P k ( T;R 3 ) and so (17.41) implies 

Wh = 0. This variable is introduced in order that (17.41) has the form of a standard 
mixed system. The discrete Helmholtz decomposition of [I] Lemma 5.4] proves 
that for the lowest order discretization k = 0, ph is a Crouzeix-Raviart function and 
so (17.41) can be seen as a generalization of the non-conforming Crouzeix-Raviart 
FEM to higher polynomial degrees. 

The inf-sup condition follows from VL^(T) C Yn fc(^T) and Curl k{7) C 
Xh(7). This and the conformity of the method lead to the best-approximation 
result 


Wp-PhWma.) + ||Curl(a-a A )|| i 2 (n) + ||V(«7 - w h )\\ L 2 {n) 

~ ( n L il3 . min ||Curl(a - Ph)\\ L 2 {Q ) 

K QhE-Xh(J) Ph€r N ,k( L) 


+ min 

s h eu h ( T) 
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Since w = Wh = 0, this is equivalent to 
lb-P/ilUa(n) + l|Curl(a - a h )\\ L 2 m 

£ ( HP “ Qh\\tf(ti) + a min ||Curl(a - Ph)\\ L 2 (n] ) • 

X( lh€X-h(J) P/i€*N,fcW 7 

The following proposition states a projection property similar to Lemma 13.121 for 
the two-dimensional case. To this end, define 

Z h {7) := {P h £ Y N , k (7) | \tv h £ U h {7) : (P h ,Vv h ) L 2 {n) = 0)}, 

W h {7) :={q h £ X h {7) \ V(3 h £ Z h {7) : (q h , Curl&) L 2 (n) = 0}. 

Since XU k (7) is the kernel of Curl : Ypj,fc(^0 —> X k (7), it holds 

CurlT N , fc (T) = Curl Z h (7). 


This implies 

W h {7) = {q h £ X h (7) | V& € Y N , k (7) : (q h , Cml(3 h ) L 2 {n) = 0}. 

Lemma 7.1 (projection property). Let q £ L 2 (Li;R 3 ) with (q, Cuvip)^^ = 0 
for all f3 £ H(Curl, LI) (that means that q is a gradient of a H((Li) function). Then 
LLj^/g-m £ Wh(7). If 7* is an admissible refinement of 7, then Hx h (7)Wh(7 i< ) C 
W h ( 7). 

Proof. Since Curl Lnj^T) C X k ( T) and Y^ tk (7) C H(Curl, LI), the assertion fol¬ 
lows with the arguments in the proof of Lemma [3.121 □ 

7.2 Adaptive algorithm 

This subsection outlines the proof of optimal convergence rates for Algorithm 16. II 
in 3D driven by the error estimators A and /U defined by the local contributions 

X 2 (7 e ,T) := WhjCurl^phWl^+hT II \Ph A ^e]e\\ 2 L 2 {e) , 

Ee£(T) 

fj, 2 (T) := ||<^ - n Xh (7)v\\ 2 L 2 { T) 

and m- Here, £(T) denotes the faces of a tetrahedron T £ 7 and hj £ Pq (T) 
denotes the piecewise constant nresh-size function defined by hj \t '■= hr := 
meas 3 (T) 1 / 3 . The refinement of triangulations in Algorithm 16.II is done by newest- 
vertex bisection m- Let T (N) denote the space of admissible triangulations with 
at most N tetrahedra more than To- As in Subsection 16.11 define the seminorm 

l(P.«^)U 3 : = sup N s ini (\\p-U Xh{7) p\\ L 2 {n) 

Ne N 0 reJL(JV) v. 

+ ll Curl ( a “ + \\<P - n x h (T)^llL2(n))- 

Assume that Assumption 16.51 holds. The following theorem states optimal conver¬ 
gence rates for Algorithm 16.11 for 3D. 
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Theorem 7.2 (optimal convergence rates of AFEM for 3D). Let s > 0. For 
0 < pb < 1 and, sufficiently small 0 < k and 0 < 9 < 1, Alaorithm \6. 1\ computes 
sequences of triangulations (T^gpj and discrete solutions (pi, o^gfsj for the right- 
hand side (p of optimal rate of convergence in the sense that 

(card(T £ ) - card(T 0 )) a (||p-p*|| L 2 (n) + ||Curl(a - a<)|| L 2 (fl) ) < \(p,a,ip)\ As . 

The proof follows as in Section [ 6 ] from (A1)-(A4) and (B) from [16] and the 
efficiency of A and p. The proof of efficiency follows with the standard bubble- 
function technique |43j . The proofs of the axioms (A1)-(A4) and (B) are outlined 
in the following. 

The axioms (Al) stability and (A2) reduction follow as in Subsection 16.21 with 
triangle inequalities, inverse inequalities, a trace inequality similar to | 10 l p. 282], 
and the mesh-size reduction property /t 3 \t < /i 3 |-t/ 2 for all T G T*\T. However, 
for (A3) quasi-orthogonality and (A4) discrete reliability, the interpolation oper¬ 
ator of } 39] cannot be applied directly to 77 * G An 7 , (7*) as done in the proof of 
Theorem 16.91 because Yn*,(T*) $7 i7 1 (D;M 3 ). This can be overcome by a quasi¬ 
interpolation based on a quasi-interpolation operator from (38] and a projection 
operator from [44] , Its properties are summarized in the following theorem. 

Theorem 7.3 (quasi-interpolation). Let T* be an admissible refinement of 7 and 
define K(T,T*) := {T G 7 \ 3Ki G 7 \ 7±3K 2 G 7 with K\ n K 2 + 0 and T n 
K 2 0}. Let 737 G Z h (f7f). Then there exists 77 G In^T), p G H 1 ^), and 
4- G AT 1 (0; M 3 ) with 


70* - 7t = Vp + 4>, 

(7T* - 7t)|t = 0 for all T G 7 \ 31(7, T*), 

II^T 1$ l|L 2 (n) + ||V$||i2(n) < ||Curl7^|| L 2( n ). 

Proof. This follows as in the proof of [441 Theorem 5.3] and with the ellipticity on 
the discrete kernel from (21 Proposition 4.6]. □ 

The differences between the proof of (A4) discrete reliability and the proof 
of Theorem 16.91 are outlined in the following. Let (py^a^) G X^( T*) x Zf l (7 ir ) 
and (pj, 07 ) G Xh(7) x Zh( 7) denote the discrete solutions to (17.31) . As in the 
proof of Theorem 16.91 let 07 ^ G Wh(7 it ) and 77 * G Zg( T*) such that p 7 — p 7t = 
<737 + Curl 77 ^. The first term of the right-hand side of 

Wpt -PT *|| 2 = (pr -Pt*,ctt*)l 2 (Q) + (po- -PT*,Curlr T J L 2 (f2) 

is estimated as in the proof of Theorem 16.91 while for the second term, the quasi- 
interpolant 77 G Yn^T) of 77 * with r<j t — 77 = Vp + 4> for p G H l (Ll) and 
4> G iL 1 (D;M 3 ) from Theorem 17.31 is employed. This yields 

(PT-PT*, Curl 77 ^ 2 ( 0 ) = (pr, Curl(r T * - r T )) i2(n) = (p 7 , Curl<f>) L2(Q) . 

A piecewise integration by parts and the arguments of the proof of Theorem 16.91 
conclude the proof. The crucial point is that 4> G iL 1 (D;M 3 ) is smooth enough to 
allow for a trace inequality. 
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(a) (b) 


Figure 1: Red-refined triangle and initial mesh for the L-shaped domain. 


The proof of (A3) quasi-orthogonality follows as in the proof of Theorem 16.111 
with the projection property of Lemma [TT] and the following modifications in (16.111) . 
Since (in the analogue notation as in (16.1111 ) a£_i e z h {7z_i) c y N>fe (T M ), 
there exists 7 m £ Z/,(Tm) with CurlyM = Curla!£_i. Theorem 17.31 guaran¬ 
tees the existence of fie-i £ T/v,fc(T^_i), p £ and <f> £ AT 1 (O;M. 3 ) with 

— 7 M — Pi-i = Vp + <L. This implies in (16.111) that 

(Curl(a M - ai-.i),pe-i) L 2 ( n) = (Curl(a M - 7m - A?-i),P 0 -i)L 2 (f 2 ) 

= (Curl$,p£_i) i2(n) . 

Since <F £ iL 1 (J7;R 3 ) is smooth enough, a piecewise integration by parts and the 
arguments of the proof of Theorem 16.91 then prove 

(Curl(a M - ae-i),pe-i) L 2 (n) < (A £ _i + pe-i) ||Curl(a M - a£_i)|| L2(n) . 

This and the arguments of Theorem 16.111 eventually prove the quasi-orthogonality. 

8 Numerical experiments 

This section presents numerical experiments for the discretization (13.4p for k = 
0,1,2. Subsections 18.1118.31 compute the discrete solutions on sequences of uni¬ 
formly red-refined triangulations (see Figure [Tal for a red-refined triangle) as well 
as on sequences of triangulations created by the adaptive algorithm 16.11 with bulk 
parameter 0 = 0.1 and k = 0.5 and p = 0.75. The convergence history plots are 
logarithmically scaled and display the error \\p — Ph\\L 2 (tt) against the number of 
degrees of freedom (ndof) of the linear system resulting from the Schur comple¬ 
ment. The underlying L-shaped domain P := (—1, l ) 2 \ ([0,1] X [—1,0]) with its 
initial triangulation is depicted in Figure [Tbl 

8.1 L-shaped domain, I 

The function u given in polar coordinates by 

u(r,<f>) = r 2 / 3 sin((2/3)^>) 
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is harmonic. For the following experiment we choose ip = 0 and uo '■= gu with 
perturbation function g G H 2 (Q), 


g(x) 


'o if \x\ < 1/2, 

< 16|x| 4 - 64|x| 3 + 88|x| 2 - 48|®| +9 if 1/2 < \x\ < 1, 
, 1 if \x\ > 1, 


such that g|r = 1 for T := dQ \ ({0} x (—1,0) U (0,1) x {0}). Since n|gQyr = 0, it 
holds UD\dti = u - Let. B 1/2(0) := {x G M 2 | |x| < 1/2} denote the ball with radius 
1/2 and midpoint (0,0). Since 5 |_b 1/2 (o) = 0 and u G H 2 (Q \ S 1 / 2 (0)), it holds 
ud G H 2 {VL). 

For non-homogeneous Dirichlet data, the jump [ph]E' r E is defined for boundary 
edges E G £, E C Te, with adjacent triangle T + by 


\Ph]E ■ t~e ■= Ph \t + ■ te — Vud ■ te ■ 

The error estimator A is then defined by (I6.2H (I6.3D . The local data error estimator 
contributions read 


p 2 (T) := ||(</> - Vu D ) - n - Vu D )\\ 2 L 2 {T y 

The global error estimator p is defined by ([6.3ft . 

The errors and error estimators for the approximation p^ G Pfc(T;R 2 ) of Vri 
for k = 0,1,2 are plotted in Figure [2] against the number of degrees of freedom. 
The errors and error estimators show an equivalent behaviour with an overestima¬ 
tion of approximately 10. Uniform refinement leads to a suboptimal convergence 
rate of /r 2 / 3 ~ ndof -1 / 3 for k = 0,1,2. The adaptive refinement reproduces the 
optimal convergence rates of ndof^ fc+1 L 2 for k = 0,1,2. Figure [3] depicts three 
meshes created by the adaptive algorithm for k = 0, 1, and 2 with approximately 
1000 degrees of freedom. The singularity at the re-entrant corner leads to a strong 
refinement towards (0,0), while the refinement for k = 0,1 also reflects the be¬ 
haviour of the right-hand side, i.e., one also observes a moderate refinement on 
the circular ring {x G \ 1/2 < |x| < 1}. The marking with respect to the 
data-approximation (p 2 > k\ 2 in Algorithm 16.ip is applied at the first 7 (resp. 5 
and 10) levels for k = 0 (resp. k = 1 and k = 2) and then at approximately every 
third level. 


8.2 L-shaped domain, II 

For / = — 1 and ud = 0 define ip(x,y) := (1/2 )(x,y) with — divy) = /. 

The error estimators are plotted against the degrees of freedom in Figure [4] for 
k = 0,1, 2. The error estimators show for fc = 0,l,2a suboptimal convergence rate 
of /r 2 / 3 ss ndof -1 / 3 for uniform refinement. The adaptive algorithm l6. ll recovers the 
optimal convergence rate of ndof~( fc+1 L 2 . Adaptively refined meshes are depicted 
in Figure [ 0 ] for approximately 1000 degrees of freedom. The strong refinement 
towards the singularity at the re-entrant corner is clearly visible. The smoothness 
of <p G Pi(fl;M 2 ) implies that the data-approximation error estimator p( vanishes 
on all triangulations for k = 1,2. For k = 0, pi does not vanish, nevertheless, 
since p 2 < k\ 2 for all £, only the Dorfler marking is applied. 
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•\A 2 + m 2 for fc = 0, uniform 
’ IIP “ Pfc II for fe = 0, uniform 

+ for fc = 0, adaptive 

Up — ph II (r2) for k = 0> ada P tive 

+ for fc = 1, uniform 

• Up - Ph II L 2 {n) for k = 1 ' uniform 

v^+ for fc = 1, adaptive 

• Up - Ph II l 2 ^) for k = 1 ’ ada P tive 

+ for fc = 2, uniform 

• Up - Phil £,2^) for fe = 2, uniform 

+ for fc = 2, adaptive 

• IIP - Ph II £2^) for fe = 2, adaptive 



ndof 


Figure 2: Errors and error estimators from Subsection 18. 11 





-1 - 0.5 0 0.5 1 


-1 - 0.5 0 0.5 1 


-1 - 0.5 0 0.5 1 


Figure 3: Adaptively refined triangulations for the experiment from Subsection 18. 11 


y* 2 

sj>? +M 2 

V ^ 2 + „ 2 

v/^ 2 +^ 2 

y^ 2 +m 2 

i/^ 2 +^ 2 


for fc = 0, uniform 
for fc = 0, adaptive 
for fc = 1, uniform 
for fc = 1, adaptive 
for fc = 2, uniform 
for fc = 2, adaptive 



ndof 


Figure 4: Error estimators for the experiment from Subsection 18.21 
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Figure 5: Adaptively refined triangulations for the experiment from Subsection 18.21 


8.3 Singular a 


This subsection is devoted to a numerical investigation of the dependence of the 
error \\p — Ph\\L 2 (n) on the regularity of a. The exact smooth solution u € C°°(Q) 
of 


—Alt = 2 sin(7rx) sin(7ry) in hi and u\r D = 0 

reads u(x,y) = sin(7rx) sin(7T7/). Define p = Vrt+Curl(5) with a £ H 1 (Fl)\H 2 ( 0) 
defined by a (r, </») = r 2 / 3 sin(2<?i/3). Then p £ H( div, 0) with — di\ p = /. 

The errors and error estimators are plotted in Figure [6] against the number of 
degrees of freedom. The convergence rate on uniform red-refined meshes for k = 
1,2 is /i 2 / 3 ~ ndof -1 / 3 and, hence, the convergence rate seems to depend on the 
regularity of a. The errors and error estimators show the same convergence rate. 
Figure [7] focuses on the results for k = 0 and uniform mesh-refinement. The error 
\\p — Ph\\L 2 (n) an d the error estimator y/ A 2 + p 2 show a convergence rate between 
h and h 2 / 3 , while ||Curl(a — u/OIIl 2 ^) converges with a rate of h 2 / 3 ~ ndof -1 / 3 
due to the singularity of a. This numerical experiment suggests that the error 
\\p — Ph\\L 2 {n) does not depend on the regularity of a (at least in a preasymptotic 
regime). The triangle inequality implies ||Curl(o; — ah)\\L 2 (n) — \\P ~Ph\\L 2 (n) + P- 
This upper bound is also plotted in Figure [7] 

Figure [8] depicts adaptively refined meshes for k = 0,1,2 with approximately 
1000 degrees of freedom. The singularity of a leads to a strong refinement towards 
the re-entrant corner. The marking with respect to the data-approximation (/i 2 > 
in Algorithm 16.111 is only applied at levels 1-5, 7, 12, and 18 for k = 0. All 
other marking steps for k = 0,1, 2 use the Dorfler marking (p 2 < kX 2 ). 
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~ y/\* + 

- Up - Phil 


for k = 0, 

p h\\ L 2 m '° rk 

A 2 -f- /a 2 for k = 0, 

IIP - PfrH I ,2 ( f J) fork 

-\J A 2 4- /a 2 for fc = 1, 

IIP - Phll L 2( n) <= 

f or k = 1 , 

IIP - Phll L 2 (n) ^ k 

^ j A 2 -f- /a 2 for k = 2, 

HP “ P hW L 2 (fl) '° ,k 

-\J A 2 + /a 2 for k = 2, 

IIP - Phil 1 , 2 ( 0 ) for k 


uniform 
= 0, uniform 

adaptive 
= 0, adaptive 

uniform 
= 1, uniform 

adaptive 
= 1, adaptive 

uniform 
= 2, uniform 

adaptive 
= 2, adaptive 



ndof 


Figure 6: Errors and error estimators for the experiment with singular a from Subsec¬ 
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Figure 7: Errors and error estimators for the experiment with singular a from Subsec¬ 
tion [531 and uniform refinement. 





Figure 8: Adaptively refined triangulations for the experiment from Subsection 18.31 
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